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We study the magnetic field dependence of the entanglement entropy in quantum phase transition 
induced by a quench of the XX, XXX and the LMG model. The entropy for a block of L spins 
with the rest follows a logarithmic scaling law where the block size L is restricted due to the 
dependence of the prefactor on the quench time. Within this restricted region the entropy undergoes 
a renormalization group (RG) flow. From the RG flow equation we have analytically determined 
the magnetic field dependence of the entropy. The anisotropy parameter dependence of the entropy 
for the XY and the LMG model has also been studied in this framework. The results are found 
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' any quench. 

O ' 

keywords: entanglement entropy, renormalization group, magnetic field; 
Pacs nos: 03.67.Bg, 03.67.Hk, 03.65.Ud, 73.43.Nq; 
f : E-mail : banasri@isical.ac.in; *: E-mail: b_ pratul@yahoo.co.in 
/: E-mail: majumdar. priyadarshi@yahoo.com 



1 



I. INTRODUCTION: 



It is well known that the entanglement entropy in quantum phase transition (QPT) of one 
dimensional spin systems in general follows a scaling law. In fact the entanglement entropy 
of a block of L spins with the rest of the system follows a logarithmic scaling law where the 
prefactor is determined by the central charge of the relevant conformal field theory. Critical 
ground states are characterized by an entropy Sl that diverges logarithmically with L having 
the relation with a coefficient given by the holomorphic and antiholomorphic central charge 
of the conformal field theory [l| . Indeed this expression corresponds to the geometric entropy 
for a conformal field theory derived by Holzhey, Larsen and Wilczek [2J. It has been shown 
in some earlier papers jj], 4] that the central charge in conformal field theory corresponds 
to the Berry phase factor acquired by a spin state when it evolves in a closed path. The 
entanglement entropy of a pure state can be reduced to the measure of entanglement of 
formation in a mixed state given by concurrence {sj. In a spin system the concurrence C 
for the entanglement of two nearest neighbor spins is found to be given by the Berry phase 
factor <f> where the Berry phase acquired by a spin state when it evolves in a closed path is 
e i<t> = e «2jr0 Jg-8|_ it m ay be noted that the central charge in conformal field theory satisfies a 
renormalization group (RG) flow as pointed out by Zamolodchikov This implies that the 
Berry phase factor </> as well as the concurrence for an entangled spin system also satisfies the 
RG flow equation. This essentially corresponds to the fact that the entanglement entropy 
undergoes a RG flow |10| . 



In some recent works 



11 



12| it has been pointed out that the entanglement entropy in 



QPT in one dimensional spin systems induced by a quench also satisfies a scaling law in a 
restricted sense such that there is a constraint on the block size L depending on the quench 
time. Indeed in this case the prefactor has a dependence on the quench time. An interesting 
result observed in this case is that the entanglement entropy in QPT in all one dimensional 
spin systems induced by a quench follows a universal behavior. In the scaling region the 
entanglement entropy undergoes a RG flow. 

It may be noted that the external magnetic field is the control parameter in QPT. In 
transverse Ising model we have a sharp critical point when the external field parameter 
A takes the value A = Ac = 1. For A > 1 the system is in a paramagnetic state and for A < 1 
the system transits to a ferromagnetic state having all the spins either in the up or down 
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direction. However for the XX model and the Heisenberg spin chain (XXX model) we have a 
critical region where the system is gapless. In fact in this case < A < 2 corresponds to the 
critical region. For A = the system attains the maximum entropy and as A increases the 
entropy decreases when at A = 2 it vanishes. In fact at A = 2 both the XX and XXX systems 
transit to ferromagnetic states. An analogous behavior appears in the Lipkin-Meshkov-Glick 
(LMG) model [k| when in the isotropic case criticality corresponds to the region < A < 1. 
Here also in the isotropic case the entropy is maximum at A = and as A increases entropy 
decreases and finally it vanishes at A = 1. For A > 1 the ground state is a fully polarized 
product state. 
In some recent works 



11 



12] we have computed the entanglement entropy in QPT induced 



by a quench for the XX, XXX and LMG models at A = 0. From an analysis of the RG 
flow we shall study here the behavior of the entropy with the increase in A until it vanishes. 
Also we shall extend our study for the XY and LMG models by quenching the system across 
quantum multicritical points by approaching along a linear path. In this case the external 
field is a linear function of the anisotropy parameter 7. We shall study the behavior of the 
entropy with the change in the anisotropy parameter. 

In secJITlwe shall consider the behavior of the entropy with the change in the external field 
(A) for the XX and XXX models. In sec lIIII we shall consider the LMG model. In sec JIVI 
we shall study the behavior of the entropy at criticality in the XY model as well as in the 
LMG model with the change of the anisotropy parameter. 



II. ENTANGLEMENT ENTROPY IN THE CRITICAL REGION OF XX AND 
XXX MODELS: 

The XX model is given by the Hamiltonian 

H = - E + + A E C 1 ) 

i i 

The criticality of the system has two-limit behavior. At A = 2 the system corresponds to 
a ferromagnetic state while at A = the system falls into the free boson universality class. 
The interval between these two points corresponds to the critical region. The entropy is 
maximum at A = and with the increase of the magnetic field the entropy decreases and 
finally at A = 2 it vanishes. At A = 2 the system transits to the ferromagnetic state when 
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the ground state corresponds to the product state. At A = the entropy scales like 

Sl ~ < —^-log 2 L (2) 

which corresponds to the entropy of a block of L spins with the rest of the system. Here c (c) 
is the holomorphic (antiholomorphic) central charge of the relevant conformal field theory 
and in the bosonic case we have c — c—1. Thus we have 

Sl ~ \log 2 L. (3) 

It has been shown in earlier papers jjj, 4| that the central charge in conformal field theory 
is related to the Berry phase factor 0, the phase being e l2n< ^ which is acquired by a spin 
state when it evolves in a closed path. It may be pointed out that the measure of entangle- 
ment given by concurrence of nearest neighbor spins in a mixed state is equivalent to__the 



entanglement entropy in a pure state [5(. It has been shown in some earlier papers js-|8] 
that the concurrence C corresponding to the entanglement of two nearest neighbor spins is 
related to the Berry phase factor and we have C = <fi. Now we note that in view of the 
relation of the central charge c with the Berry phase factor <p Zamolodchikov's c-theorem 9] 
representing the RG flow of the central charge can be transcribed in terms of 6 and the 

n 

concurrence C which essentially implies the RG flow of the entanglement entropy [10[ . The 
RG flow suggests that the entropy decreases along the flow and we have 

4,0, (4) 
where L is a length scale. From this we have 

w a In L = alog 2 L. (5) 

L 

For a pair of nearest neighbor spins (L = 2) we have a = <p which corresponds to the 
concurrence of two nearest neighbor spins in the system. So for the entanglement entropy 
for a block of L spins with the rest we write 



S, 



log 2 L. (6) 



At the critical point A = 0, the system belongs to the boson universality class and the Berry 
phase factor which is identical with the concurrence for a pair of nearest neighbor spins 
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in an antiferromagnetic system is given by C = <p = 0.386 1 151 ] . This is very close to 
the prefactor 1/3 in ([3]) derived from the conformal field theory. It may be mentioned here 
that the entanglement of a block of L spins with the rest of the system can be considered 
to be equivalent to the entanglement between a single spin representing the block spin with 
another spin represented by the rest of the system in block variable RG scheme. In view 
of this Sl in (EJ) can be considered as the concurrence C for the entanglement between the 
pair of this two block variable renormalized spins in a mixed state. The slight departure of 
the prefactor 1/3 in ([3]) from the value = 0.386 in may be associated with the block 
variable renormalization of the spi n sytem which induces change in the coupling constant. 
In fact in our earlier works [ll[ Q] from an analysis of the transverse Ising model we have 
introduced a correction factor 0.926 associated with the block spin variable. 
The introduction of a quench incorporates a new length scale given by the Kibble-Zurek 
(KZ) correlation length £ If]- 20], which scales like £ ~ y/r, r being the quench time 12 ]. 



Taking into account this aspect we have founc 
spins with the rest at A = for the XX model 



the entanglement entropy for a block of L 
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S L (A = 0) 



l °9i L loq 2 L 

x 0.926 « 3.7—^-. 

log 2 i l °92r 



(7) 



It should be mentioned that the value of L is here restricted and therefore the entropy does 
not rise with the addition of spin in the block size indefinitely. In fact from the constraint 
Sl{t)/ S max where the maximum value of the entropy is given by [12 ] 



Smax = 2(\(f)\log 2 l + 1) x 0.926 « 0.25 Inr + 1.85, 



(8) 



we have the relation 



In L < 0.07(Zn r) 2 + 0.5/n r. 



(9) 



Within this restricted region the entropy undergoes the RG flow. 

Now to study the magnetic field dependence of the entropy in the critical region < A < 2 
the time dependent magnetic field is taken to be given by 



2t 

X(t <0) = 2-—. 

T 



(10) 



so that at t — r the system resides at the critical point A = and it evolves toward A = 2 
when at the end t — it reaches there. Now transcribing the RG flow equation in terms 



of time (L = ct), we write 

fi " <0, (11) 



which implies 

a In t/r = alogi t/r. (12) 



t 



we have 

^log 2 t/r. (13) 













<P 




t 





As before identifying a 



However there is a caveat here. When we consider that the time dependent magnetic field 
X(t < 0) traverses the critical region < A < 2 starting from A = at t = r and ending at 
A = 2att = 0ina closed circuit it is noted that a spin state while traversing this closed 
path will acquire the Berry phase factor \<p\. However when we transcribe the RG eqn.(TjI|) 
in terms of t/r as shown in (llip we note that t/r is restricted in the region < t/r < 1. 
So when we map the closed circuit < A < 2 onto the circuit < t/r < 1 the Berry phase 
factor acquired by a spin state while traversing this closed path will be half of that acquired 
in the former case. So for the effective phase factor we write 

\4>\eff = \\^ (14) 

and we have 

AS = |0| e// log 2 t/r, 

= 1$ log 2 t/r. (15) 

Now from ( FlUj) we have for the entropy variation with A from the value at A = 

AS(A) = i|#<^(l-0. (16) 

As mentioned above the Berry phase factor <fi which corresponds to concurrence for a pair 
of nearest neighbor spins in the spin chain which at A = corresponds to that in the 
antiferromagnetic system is given by <fi = 0.386. Incorporating the correction factor 0.926 
we have the prefactor 0.35 which is very close to the value 1/3 derived from conformal field 
theory. So we write 

S L (X) = S L {\ = 0) + ^log a (l-± 

« S L (X = 0) + ^logJl-^\. (17) 



This is valid for all values of L which satisfy the constraint (Q. It may be noted that an 
analytical expression for the magnetic field dependence of the entanglement entropy has 
been given by Jin and Korepin 



2l|. It has been observed that their result is compatible 
with the numerical studies which fixes the value of the constant term to be added to the 
expression (j 



22|. 



It is observed that the A dependent term in expression (TTTjl is independent of r and thus will 
be valid for QPT without a quench. Indeed it is found to be in excellent agreement with the 
numerical results obtained by Lattore and Riera as shown in fig. 2 in 22]. For QPT induced 



by a quench taking the expression for Sl (A = 0) as given by ([7]) we can compute Si, (A). 
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FIG. 1: (Color online) Figure (left) shows S L vs. A for L = 50 and r = 200, 400, 800. Figure (right) 



depicts the variation of Sl with L for r = 200 anc 
results for A = and A = 1.9 with that derived in 



A = 0, 0.5, 1.0, 1.5, 1.9. We have compared our 
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In fig{H (left) we plot the variation of Sl vs. A for a fixed value of L = 50 and for different 
values of r = 200, 400, 800. In figJU (right) we plot S L vs. L for a fixed value of r = 200 
and for different values of A = 0, 0.5, 1.0, 1.5, 1.9. We have compared our results with that 
obtained by Lattore and Riera 22| for A = and A = 1.9 without any quench in the 
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thermodynamic limit. From the figiU (right) we note that for r = 200 the results are found 
to be in very good agreement with their findings without the introduction of any quench. 
For higher values of r the entropy decreases as follows from (j7j). 
The Hamiltonian for the XXX model is given by 

H XXX = K^f+i + °Mu + + A £ <■ (18) 

i i 

The critical behavior of the model is analogous to that of the XX model. It has two-limit 
behavior. At A = 2 the system represents the ferromagnetic state and at A = the system 
corresponds to the antiferromagnetic state. The interval < A < 2 is gapless and hence 
critical. As at A = this corresponds to the antiferromagnetic state the entanglement 
entropy for a block of L spins with the rest of the system is identical with that of the XX 
model. This is valid for the quench induced QPT also. Just as the XX model here also the 
entanglement entropy is maximum at A = and decreases with the increase in A until at 
A = 2 it vanishes. The variation of Sl(A) with different values of A will be identical with 
that of the XX model. 

III. ENTANGLEMENT ENTROPY IN THE CRITICAL REGION OF THE 
LIPKIN-MESHKOV-GLICK MODEL: 

The Hamiltonian for the Lipkin-Meshkov-Glick (LMG) model is given by 

H = ^ E + -rW) + A E ( 19 ) 

i<j i 

N being the total number of spins. In contrast to the conventional spin model in the LMG 
model each spin interacts with all the spins of the system with same coupling constant. This 
introduces the loss of the notion of geometry as there is no distance between the spins. Thus 
we cannot consider the notion of a block of spins as a set of contiguous spins here. The 
symmetry of the Hamiltonian suggests that the ground state belongs to a symmetric subspace 
where all the spins are indistinguishable and this subspace restricts the entanglement entropy 
of a block of L spins with the remaining spins. However the scaling behavior shows a similar 
pattern with that of the XX model where conformal symmetry plays a significant role in the 
critical region. However in the LMG model though the scaling law is similar to that of the 
XX model it has nothing to do with the underlying conformal symmetry. 
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In the isotropic case with 7 = 1 the Hamiltonian can be written in terms of the total spin 
operator S a = 1/2 as 



H 



(S z ) 2 - — ) +2XS Z . 



N 



(20) 



Now as shown in an earlier paper [12j if we introduce point-splitting regular izat ion so that 
we write 

S 2 = S k .S' k S kk i, (21) 

where k and k! are two adjacent sites with an infinitesimal distance k — k! = e and the 
relation ( 121 j) is satisfied in the limit e — > 0, we can take 



na \ ^ a riot _ ^ \ ^ q 



(22) 



2^ *' " 2 

with i and j being two adjacent sites with an infinitesimal distance. Considering only nearest 
neighbor interactions we can take the Hamiltonian in the regularized form 12| 



H 



reg 



1 

2iV 



L hi 



2N \ ^ 1 1 ^ 

\ hi i 



07 



1. 



(23) 



From this we note that when we consider 1/2N as a dimensionless coupling constant J we 



can write H, 



reg — H 1 + H 2 - 1 with 
Hi = J 



v r ( ^x ^x 1 _i/ _y 1 _*_je\ , ^ \ r _2 



A 

2J 



0", 



(24) 
(25) 



hi, = -./ ! \ a: a] - — X). 

*j i 

The Hamiltonian given by ( 1231) can be split into various other H\ and i?2 such that both of 
them describe critical systems. However as the LMG model is characterized by the fact that 
each spin interacts with every other spin so that for even (odd) number of interacting spins 
we have bosonic (fermionic) systems, we have chosen the regularized Hamiltonian such that 
one of them represent bosonic and the other fermionic features. From (j24p and (I2"5"j) we note 
that Hi and H 2 represent the XXX model and Ising model respectively. 
To consider the criticality of the system we observe from ( 1231) . ( 1241) and (1251) that this is a 
combination of the XXX model given by Hi and the Ising model with magnetic field along 
the negative z direction given by H 2 . It is observed that in this regularized form A/2J 



denotes the intensity of the magnetic field. The XXX model is characterized by the fact 
that for 2 > |A|/2J > the system is gapless and hence critical. For |A|/2J = 2 we have 
|A| = 2/N as we have J = 1/2N. Since for an entangled spin state the minimum number of 
spins must be 2 we find that at criticality |A| lies in the interval < |A| < 1. For the Ising 
chain given by H 2 the region |A| > 1 corresponds to the fact that all spins are polarized 
along the negative z-axis. In the interval < |A| < 1 as |A| is tuned from 1 to the spin 
system will undergo a transition when down spins will be excited so that at A = all spins 
will settle down with opposite orientation. In between these two points with the tuning 
of A spins evolve through a paramagnetic state. Thus during the transition in the interval 
< |A| < 1 spins evolve through a situation which is similar to that of the transverse 
Ising model. It is observed that when the LMG model is recast in the regularized form the 
fact that the critical region manifests the same logarithmic scaling law as observed in XXX 
model and the transverse Ising model can be understood from the conformal symmetry in 
the critical region in these systems. Indeed the point-splitting regularization unveils the 
underlying conformal symmetry at criticality in this system which is lost in the sharp point 
limit. 

In the regularized Hamiltonian introducing a quench we consider that the system transits 
from the point A = where the entropy is maximum towards the value |A| = 1 where the 
entropy vanishes. So for QPT induced by a quench we consider that the time dependent 
magnetic field behaves as 

A(t < 0) = 1 - - (26) 

T 

so that at t = r, r being the quench time, A = and at the final state, when t = we 
have A = 1. Now from the RG flow (j4]) and transcribing it in terms of t/r we have for the 
variation of the entropy with A from the value at A = 

AS « \4>\log 2 t/r, (27) 



with cf) = \<f)\xxx + 1 4>\ ising- It has been shown in an earlier paper [TjJ that the Berry 
phase factor which is equivalent to the concurrence for the nearest neighbor spins in the 
transverse Ising model is given by C = \4>\i s ing — 0.18. For the XXX model we have 



\cf)\ xxx = 0.386 [14], q4- Thus \<f>\ = 0.386 + 0.18 = 0.566. Taking into account the 
correction factor 0.926 as mentioned in the previous section we have the effective value of 
\4>\ e ff = 0.52. Now from fl26|) we have t/r = 1 — A. However we have to take care of the fact 
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that in the regularized Hamiltonian for the XXX model the external field is in the positive 
z-direction while for the Ising model this is in the negative z-direction. So in the sharp point 
limit we have to take both these orientations of the magnetic field. Thus from (j27p we write 
in the sharp point limit for the entropy variation with A from the value at A = 

AS(X) = $\ eff [log 2 (l-X)+log 2 (l + X)}, 

= 0.52 log 2 (l - A 2 ). (28) 

Thus we have 

S L (Xn = 1) = S L (X = 0,7 = 1) + 0.52 log 2 {\ - A 2 ), (29) 



which is in excellent agreement with the value obtained numerically 23 1 

S L (A, 7 = 1) = S L (X = 0,7 = 1) + \ log 2 (l - A 2 ). (30) 

It may be mentioned here that when QPT is induced by a quench the value of the en- 
tanglement entropy for the isotropic LMG model (7 = 1) at A = is found to be given 



by 



12| 



^>l{t) ps 2— - x 0.926 3.7- . (31) 

<j>log 2 £ lo 92r 



From the constraint Sl(t)/ S max < 1 we find the constraint for L [Tj 

In L\ LMG < 0.097(/n r) 2 + 0.5 In r. (32) 

Within this restricted region of L the entropy scales like log 2 L. Thus within this restricted 
region of L the RG flow will be valid. So from (1281) we can write explicitly 

S L (X, r, 7 = 1) = S L (X = 0, r, 7 = 1) + 0.52 log 2 {\ - A 2 ), (33) 

with 

&(A = 0, T , 7 =1>»^, (34, 



where L is constrained by the relation ( 1321) . 

In figj2] (left) we plot the variation of Sl with A for a fixed L = 50 and for different values of 
t = 200, 400, 800. In figj2] (right) we show the variation of Sl with L for a fixed r = 200 and 
for different values of A = 0.25, 0.5, 0.75, 0.9. The results are compared with that obtained 
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23| without introducing the quench extrapolating it to the thermodynamic limit. It is 
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FIG. 2: (Color online) Figure (left) depicts the variation Sl(A,7 = 1) with A for a given L = 50 
and for different values of r = 200, 400, 800. In fig|2] (right) we have plotted the entropy vs. L for a 
fixed value of r = 200 and for different values of A = 0.25, 0.5, 0.75, 0.9. The results are compared 
with that obtained in 23] extrapolating it to the thermodynamic limit. 



noted that with increase of r, the entropy decreases. 

In the anisotropic case (7 7^ 1) for the LMG model in the thermodynamic limit the ground 
state for A > 1 represents the product state. For A — > the entanglement entropy saturates 
and goes to a constant that depends on 7. For 7 = the ground state is degenerate and 
lives in the subspace of ±x direction. These two different phases suggest the existence of a 
QPT between A > 1 and A < 1. We now note that 7 = is the representative of the class 
7 7^ 1 and for 7 = 0, A > the regularized Hamiltonian corresponds to the transverse Ising 
model. Indeed from the relations (JT9l) and (1231) we note that for 7 = there is an extra term 
—J^2 a i a j m t ne second term on the r.h.s. of (|23|) . This modifies the Hamiltonian H 2 in 
fl25|) which is now given by the equation 



H 2 = # 2 ( 7 = 0) = -J [J2 



a i a 3 



A 



afo* - —a? 

1 J 



2J 



(35) 
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The total regularized Hamiltonian 

H = H(j = 0) = H 1 + B 2 (36) 

effectively corresponds to that of the transverse Ising model. Now from the relation (T2"6"j) we 
note that at t = it reaches the critical point A = 1 starting from A = at t = r and we 
find the entropy around criticality in the thermodynamic limit 

Sk(A,7^1)«ifofc(l-A). (37) 

The prefactor 1/6 corresponds to the prefactor associated with the entanglement entropy 
of the transverse Ising model [ll[. This is identical with the result obtained by Lattore 



et.al. 



23] from numerical studies. It is to be mentioned that just like the isotropic case for 



QPT induced by a quench, L is restricted by the constraint given by ( 1321) . 



IV. VARIATION OF ENTROPY WITH THE ANISOTROPY PARAMETER: 

We extend our study here when QPT is subject to a quench across quantum multicritical 



points by approaching along a 
of the external magnetic field 



inear path formulating the anisotropy parameter dependence 



24|. We take 



A( T ) = 1- T (t), t<0 ( 7 ^1). (38) 



The XY model Hamiltonian is given by 



H xy = -YA ^vM+i + ^ofJLx + K ) • (39) 



i + T.** 1-7 



i 

For 7 = 1 it corresponds to the transverse Ising model. For 7 ^ the system falls into the 
free fermion universality class and is critical at A = 1. For 7 = the system reduces to the 
XX system which at A = corresponds to the free boson class. To study the 7 dependence 
of the entropy we take into account the expression (138D exhibiting the 7 dependence of A 
implying multicritical behavior. When we introduce a quench we take that at t = r, r being 
the quench time the critical point 7 = 1 is reached so that we have 

7(*) = --, *<0. (40) 

T 

From (!38l) . this implies 

X(t < 0) = 1 - -, (41) 
r 
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so that at t = r we have A = 0. The critical point A = 1 is reached for the fermion 
universality class at t = by evolving from the base point A = at t = r. At the critical 
point the entanglement entropy for the transverse Ising model follows the scaling law [l| 

S L ~ hog 2 L, (42) 

which follows from the relevant conformal field theory. From the relationship between the 
Berry phase factor obtained by a spin state when it evolves in a closed path with the 
concurrence for a pair of nearest neighbor spins the prefactor is found to be 0.18 7] which is 
very close to the factor 1/6. Indeed this determines the correction factor (l/6)/0.18 = 0.926 
introduced earlier which incorporates the correction for block variable RG scheme. Now as 
discussed in secJTT]from the RG flow equation (jlj) and transcribing it in terms of t/r we have 
from (|40p the variation of the entropy S with 7 from the value of the entropy at A = 1, 7 = 1 

AS( 7 ) = ~% 2 t/r, 

= ^log 2 7- (43) 

Thus we write 

S L (X = 1, 7) = S L (X = 1, 7 = 1) + hog 2l , (44) 

6 

which is satisfied for the fermion universality class (7 7^ 0) in the thermodynamic limit. This 
result is identical with that obtained by Vidal, Lattore, Rico and Kitaev 1] from a detailed 
computational study. It is observed that the singular behavior at 7 = is due to the fact 
that at this point the system represents the XX model which at A = corresponds to the 
free boson class. 

It has been shown in an earlier paper 



ll| that when QPT is induced by a quench the scaling 
law of the entanglement entropy of the transverse Ising model (7 = 1) at the critical region 
is given by 

, . 3.7log 2 L 
Sl(t) « - 45 
log 2 r 

and the validity of this scaling property is restricted to the region given by the condition 
Sl(j) / S max < 1 which implies 

In L < 0m(ln r) 2 + 0.5 In r. (46) 
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Within this restricted region of L the RG flow equation is valid and from (1441) we can write 
for 7 7^ 

3 71nn n 1 . 1 

(47) 



3.7log 2 L 1 



/ooor 6 

n 

To compare it with the exact value of the entropy obtained numerically by Vidal et.al. [1] 
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FIG. 3: (Color online) Fig. 3 (left) shows the variation Sl with 7 for a fixed value of L = 30 and 
for different values of r = 1600, 3200, 6400, 12800. Figj3] (right) shows the variation of S L with L 
for a fixed r = 12800 and for different values of 7 = 0.25,0.5,0.75,1. The results are compared 
with that obtained in [l| for 7 = 1. 

at 7 = 1 in QPT without quench it may be mentioned that the studies in [l] involve infinite 
Ising chain. However the relation (j46j) implies that the maximum value of L is restricted 
by the quench time r and indicates that for small r, L max is very small. In view of this 
we take large r for comparison. In figj3] (left) we plot the entanglement entropy vs. 7 
for a fixed L = 30 and different values of r = 1600, 3200, 6400, 12800. In figfj (right) we 
show the entanglement entropy vs. L for a fixed r = 12800 and for different values of 
7 = 0.25,0.5,0.75, 1. It is observed that for 7 = 1 our results are in good agreement with 
that of Vidal et. al. [I] derived from numerical studies without introducing any quench. As 
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in other spin systems we note that the entropy decreases with the increase in r. 
It may be mentioned here that a specific property of QPT in transverse Ising model is 
that the maximum of the concurrence for the entanglement of a pair of nearest neighbor 
spins does not occur at the critical point A = 1 but away from it. Indeed at the critical 
point the correlation length diverges and each site develops entanglement with its nearest 
neighbor site. In view of this the entanglement of a certain pair of nearest neighbor sites is 
distributed in such a way that it saturates the constraints of entanglement sharing. Osborne 
and Nielsen |s| have conjectured that the ground state at the critical point actually saturates 
the bounds of entanglement sharing so that it is maximally entangled in this sense. This 
would imply that when the system reaches criticality the entanglement is distributed to 
more remote pairs and entanglement sharing would have to occur at the expense of the two 
party entanglement. In our present analysis we note that near criticality for the XY model 
with 7 7^ the variation of the entanglement entropy with A is given by 

S = ±\log 2 (l-\)\ (48) 

which follows from (14TI) and (|4~3|) . From this we observe that near criticality the concurrence 
for a pair of nearest neighbor spins satisfies the relation 

C = ±\log 2 (l - \)\ (49) 



Here the prefactor is taken to be half of that in ff48|) due to the fact that a spin can form a pair 
with its nearest neighbor in both left and right sides and the effect of this is incorporated in 
the expression of the entanglement entropy 3, ll2| . Thus we note that we have a logarithmic 
divergence at the critical point A = 1. This indicates that at this point the correlation length 
diverges and entanglement is shared among the remote sites. As the concurrence decreases 
along the RG flow the maximum will occur at a point A < 1 close to the critical point A = 1. 
In a similar way we can analyze the LMG model given by the Hamiltonian ( TT9|) . For 
A = 1, 7 7^ 1 we can derive the variation of the entropy with 7 from the value at A = 1, 7 = 
using the fact that 7 = represents the class 7 ^ 1 and A = 1 is the critical point. For the 
evolution of 7 with time we take the relation 7 = 1 — t/r which implies that at t = r we 
have the base point 7 = and it evolves with time so that at the end t = it reaches the 
point 7=1. From the RG flow equation (j4j) the deviation of entropy for 7 7^ 1 from the 
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value at A = 1, 7 = is found to be 



AS(A = 1, 7 ) = ^log 2 t/r, 
6 

= ^2(1-7), 



(50) 



so that we have 



S L (X = 1,7) = S L (\ = 1,7 = 0) + -log 2 (l -7). 



(51) 



This is valid for all 7 in the region — 1 < 7 < 1 . The prefactor 1/6 follows from the fact that 
at A = 1, 7 = the regularized Hamiltonian corresponds to the critical transverse Ising 
model where we have the prefactor 1/6 for the scaling law. This is consistent with the result 
obtained by Lattore et. al. [23| from numerical studies. Now to compute the entanglement 
entropy for the LMG model at 7 = in QPT induced by a quench, we note that the prefactor 
depends on the quench time which arises from the fact that the final state represents a kink- 
antikink chain with lattice spacing approximately given by the Kibble-Zurek correlation 
length £. However in the LMG model a correlation length characterizing the typical distance 
between defects cannot be introduced though we can estimate the fraction of flipped spins 
after the quench. Now from the regularization scheme for 7 = we note that the regularized 
Hamiltonians for 7 = 1 and H for 7 = differ by a term — J £V . of crj as is evident from ( 125]) 
and (j35p . This essentially introduces for the contribution of H 2 in ( 135]) an extra component 
corresponding to the transverse Ising model. Evidently the value of |0| e // i n the prefactor 
\<p\log2^ associated with the quench time in the denominator of (|3T|) will be given by 



where 0.926 is the correction factor as introduced above. This is essentially identical with 
the corresponding expression for LMG model with 7 = 1 as in ( I3TT) . This term arises from 
the fraction of defects formed during critical slowing down and the regularization scheme 
suggests that this remains unaltered irrespective of the value of 7. However as for 7 = 0, the 
total regularized Hamiltonian essentially corresponds to that of the transverse Ising model 
the scaling law of the entanglement entropy at A = 1, 7 = will be identical with that 
model viz. Ij^log^L. This implies that the numerator in ( )3T|) will be modified now and is 
given by l/6log2L. Thus for QPT induced by a quench, the entanglement entropy for the 



\Meff = \<f>\xxx + \<f>\isin 9 = (0.358 + 0.18) x 0.926 = 0.52. 



(52) 
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LMG model at 7 = can be expressed as 

£>£ (A = 1,7 = 0)\lmg 



2 ¥°92L 

0.52 log 2 ( 
37 ¥°9z L 



x 0.926 



0.52 log 2 r 



This suggests that we have the relation 



S L (A = l, 7 ^l)| iMG ~3.7- * l092L 



+ ^og 2 (l -7) 



(53) 



(54) 



0.52 log 2 r 6 

as follows from ( 15"TT) . It is to be mentioned that L is here restricted by the constraint ( 1321) . 
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FIG. 4: (Color online) Fig. 4 (left) shows the variation of the entanglement entropy with 7 for 
L = 125 and r = 1600, 3200, 6400, 12800. FigfJ (right) shows the variation of the entropy with L 
for r = 12800 and 7 = 0, 0.25, 0.5, 0.75. In figH] (left) we have compared our results with that of 
Lattore et. al. derived from numerical studies without introducing quench extrapolating it to 
the thermodynamic limit. 

In figJH (left) we show the variation of the entanglement entropy with 7 for a fixed value 
of L = 125 and different values of r = 1600, 3200, 6400, 12800. In figH (right) we plot 
the variation of the entropy with L for a fixed value of r = 12800 and different values 
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of 7 = 0, 0.25, 0.5, 0.75. In figlU (left) we have compared our results for r = 12800 with 



that of Lattore et.al. 23( and is found to be in good agreement with those results derived 
numerically without introducing quench. As it has been pointed out that for 7 ^ 1 the 
model represents the class for 7 = and for A 7^ the regularized Hamiltonian corresponds 
to that of the transverse Ising model, this explains the fact that the LMG model for 7 7^ 1 as 
A varies away from the critical value is analogous to the transverse Ising model as observed 



by Lattore et. al. 



23). 



V. DISCUSSION: 

It may be mentioned that the logarithmic scaling law was derived for geometric entropy for 
a conformal field theory [2] and this expression has been confirmed for several critical spin 
chains [l| • The present analysis suggests that the logarithmic scaling law of the entanglement 
entropy of ID spin systems is a consequence of the RG flow which follows from (J4j) - (EI) • 
However there are situations which do not obey the logarithmic scaling law. Indeed Calabrese 
et.al. 
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261 ] have pointed out that after a global sudden quench the entanglement entropy 
increases linearly with time. These authors have shown that when the system is prepared in 
a pure state \ipo) which corresponds to an eigenstate of the Hamiltonian H(\ ) with Ao 7^ A 
and at time t = the parameter is suddenly quenched from Ao to A, then the entanglement 
entropy of ID spin systems in an interval of length I increases linearly with time up to a 
certain point after which it saturates at a value proportional to / with a coefficient depending 
on the initial state. The behavior has been interpreted as a consequence of causality. In this 
context it may be added that for a quench induced QPT the logarithmic scaling behavior 
persists, though in a restricted sense, where there is a critical slowing down and the control 
parameter A is a linear function of time. In fact within the scaling region, we have the RG 
flow of the entropy. However for a sudden quench the RG equation will not be satisfied and 



hence we have a deviation 
In some earlier papers [ill . 



rom the logarithmic scaling law. 

I2I it has been pointed out that when QPT in one dimensional 



spin system is induced by a quench the entanglement entropy satisfies a scaling law with 
a prefactor which depends on the quench time. However in this case for the entropy of a 
block of L spins with the rest of the system the block size L is restricted by a constraint. 
From the RG flow equation we have derived here the variation of the entropy with external 
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magnetic field. The results are found to be in excellent agreement with that obtained from 
numerical studies by other authors without the introduction of the quench. For the LMG 
model we have taken into account the point splitting regularization and we have deduced our 
result from the regularized Hamiltonian. This unveils the underlying conformal symmetry 
at criticality which is lost at the sharp point limit. 

For the XY model and the LMG model we have studied the behavior of the entanglement 
entropy with the variation of the anisotropy parameter 7 when QPT is subject to a quench 
across quantum multicritical points by approaching along a linear path formulating the 
anisotropy parameter dependence of the external magnetic field. Using the RG flow equation 
we have considered the anisotropy parameter dependence of the entropy and is found to be 
in good agreement with that obtained by other authors flil^. For the LMG model as 7 = 
is representative of the class 7 ^ 1 and for 7 = the regularized Hamiltonian corresponds 
to that of the transverse Ising model, the surprising behavior that for 7 ^ 1 and as the 
magnetic field parameter A departs from the critical value the system is analogous to the 
transverse Ising model is well explained. 

Finally we point out that the RG flow equation for the central charge as well as for the Berry 
phase factor which is attained by a spin state when it evolves in a closed path essentially 
determines the RG flow equation of the entanglement entropy in QPT. This appears to be 
a strong tool to study the magnetic field dependence as well as the anisotropy parameter 
dependence of the entanglement entropy. 
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